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Recently there has been interest in the idea of quantum computing without control of the 
physical interactions between component qubits. This is highly appealing since the 'switching' 
of such interactions is a principal difficulty in creating real devices. It has been established 
that one can employ 'always on' interactions in a one-dimensional Heisenberg chain, provided 
that one can tune the Zeeman energies of the individual (pseudo-)spins. It is important to 
generalize this scheme to higher dimensional networks, since a real device would probably be 
of that kind. Such generalisations have been proposed, but only at the severe cost that the 
efficiency of qubit storage must fall. Here we propose the use of multi-qubit gates within 
such higher-dimensional arrays, finding a novel three-qubit gate that can in fact increase the 
efficiency beyond the linear model. Thus we are able to propose higher dimensional networks 
that can constitute a better embodiment of the 'always on' concept - a substantial step toward 
bringing this novel concept to full fruition. 

Quantum computation (QC) holds the promise of manipulating information in fundamentally new ways, 
performing tasks that are effectively impossible with any classical computer. However, the challenge of 
experimentally realizing such a device is extremely difficult. The archetypal theoretical model for QC 
assumes that one can freely perform various kinds of manipulation on the elementary qubits. (Here we 
will use the term 'qubit' to refer to the units of quantum information, and the term 'spin' to refer to the 
physical two-state systems that represent the qubits - an actual electron spin being one obvious example.) 
The required manipulations include all possible transformations on each qubit in isolation, the so-called 
one-qubit gates, together with the ability to switch 'on' and 'off' the interactions between spins in order to 
perform an entangling two-qubit gate. A literal embodiment of this model would therefore call for at least 
two different forms of control, which together can achieve arbitrary local manipulations of this kind"'^. There 
do exist many such schemes, for example the Kane proposal^ where a combination of global EM pulses and 
local eletrostatic gates is required. However, there have been several theoretical efforts to derive alternatives 
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to the archetypal model, which would be less difficult to implement experimentally. Considerable attention 
has been focused on the idea of 'all Heisenberg' computational^, where one only needs to manipulate the 
interactions between spins (dispensing with the need to separately manipulate isolated spins), at the cost of 
encoding each qubit over two or three spins. This idea is clearly a powerful simplification, however in many 
systems it is in fact the control of the spin-spin interactions that is the most difficult aspect. A related idea 
is that of global control: here one dispenses with the ability to localize manipulations and instead applies all 
manipulations globally across the entire device^. Again, this idea still requires one to switch interactions, 
albeit collectively (unless they are of the relatively passive Ising typei). 

Recently however, ideas have been advanced for performing quantum computation in strongly coupled 
systems without the need to switch interactionsSiSiifliii. It has been shown that one can perform all the 
necessary one-qubit and two-qubit gates simply by 'tuning' the Zeeman splitting of the individual qubits. 
This approach involves using some spins to hold qubits, while others are employed to act as passive barriers 
between qubits. Interestingly, this 'barrier-spin' concept proves to be compatible with the concept of global 
control, so that in one variant of the idea it is not necessary to 'target' the Zeeman tuning onto single 
spins within the device^. There is a cost associated with all existing variants of the scheme: at least one 
barrier spin is required for each qubit-bearing spin, and therefore the ratio R of qubits to spins is limited 
by i? = 1/2. This ratio is one important measure of the efficiency of an architecture, and frequently the 
QC schemes that aim to lighten experimental requirements come with the cost that R must be significantly 
less than unity. For example, the 'all Heisenberg' scheme mentioned above requires i? = 1/3 (assuming 
an isotropic environment). However, the barrier scheme has so far only been fully explored for the case of 
a one-dimensional system. This is a significant limitation: although in principle quantum algorithms can 
be performed on one-dimensional architectures without losing their advantage over classical machines, in 
practice two- or three-dimensional structures will prove considerably more effective-'^^ . There has been a brief 
discussion of the generalization of the ID results to higher dimensional networks^"'^, however the approach 
taken there comes at the cost that the proportion of inactive barrier spins must increase and consequently 
R falls to i? = 2/5 at best. It would therefore be very desirable to find a new approach to implementing 2D 
and 3D barrier networks, whereby the efficiency can at least match the ID case, if not exceed it. This is the 
problem we address in the present paper. 

We will introduce the necessary concepts by referring to the earlier work on one-dimensional arrays, 
and then continue to examine the multi-qubit nodes that are possible in higher dimensional networks, and 
which (we will show) can offer superior values of R. We will present the analysis through an explicit 
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matrix diagonalisation method that is, we hope, straightforward to comprehend for the interested non- 
specialist. In our models we will assume a general anisotropic Heisenberg interaction of the form Jzcrf cr^ + 
Jxy{<^^ '^'^X '^X)^ which includes the well known isotropic form, and the purely planar "XY" interaction. 



as special cases. There is a considerable range of pro mising physical systems associated with thi s spect rum 
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of interactions (for the isotropic limit, see e.g. Refs. 




13ll4j . and for the anisotropic case, Refs. 



Note that the XY limit is also relevant to Forster-Dexter processes, e.g. in the context of excitonic exchange 
in biological molecules. 

The Heisenberg-to-Ising Transition and the Barrier Scheme. 

Consider a chain of spins coupled via the Heisenberg interaction. Suppose that each coupled pair is 
far off-resonance, in the sense that the difference between their Zeeman energies is much greater than the 
coupling strength Jxy ■ In this limit a native Heisenberg interaction Jzcrf'^X ^ Jxy{<^i <^2 '^^X '^X) tends to 
an effective Ising coupling Jzcr^ cr^ (for a discussion of this phenomenon, see Ref. 11). This can be exploited 
as the basic mechanism for quantum gate operations, if one assumes that the Zeeman discrepancy can be 
dynamically tunec?: one would tune a set of spins into resonance to perform the gate, and then restore the 
discrepancy to return the device to its passive state. Analysis of the process is simplified if we assume very 
abrupt tuning, but this is not a necessary condition^^. In Ref. we showed how to perform this process on 
a one-dimensional chain. We introduced the idea of using a 'barrier' spin between the qubit-bearing spins, 
such that the barrier becomes entangled with the qubits only during gate operations (see Fig. 1, upper 
part). This barrier spin is initially in a known eigenstate, | t) say, and the duration of the on- resonance 
phase must be such that the barrier returns to this eigenstate. This condition corresponds to achieving 
coincidence in the revival times of different computational basis states, as will be seen presently. Previous 
work has only shown that this is possible for barriers separating two qubits. This is of course the only case 
of interest for a onc-dimensional chain (given nearest-neighbor interactions), however it is not the natural 
arrangement for two- and three-dimensional networks. 

In principle one can generate a barrier-based network by taking any arrangement of qubit-bearing spins, 
and simply introducing a barrier-spin between each interacting pair. However, a moment's reflection reveals 
that this inevitably leads to a rather large proportion of 'wasted' barrier-spins. Recall R is the ratio of 
qubit-bearing spins to the total number of spins in the device. Clearly, one would wish R to be as high as 
possible, given the difficulty of protecting spins from decoherence etc. Let us assume here that each qubit 
is represented by the state of a single spin (this docs not preclude encoding of 'logical' qubits over several 
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FIG. 1: Possible geometries for the 'barrier' QC scheme, (a) The estabhshed one-dimensional 
scheme. Letters A, B, 6 denote Zeeman energies; A and B are fixed but 6 is abruptly changed 
as shown by the graph to the right. Essentially the barrier is lowered (brought into resonance) to 
allow the neighboring qubits to interact, and is subsequently raised later when the barrier 'revives' 
into a product state. The lower illustrations show the logical extension to three-qubit (b) and 
four-qubit (c) nodes, which can be supported in higher dimensional arrays. Successive illustrations 
(vertically) correspond to the three rows in the ID case above. Fig. 2 shows how these nodes might 
be combined to form networks. 

of these elementary qubits, for the purpose of fault tolerance etc). Then for a one-dimensional array, every 
other spin is acting as a barrier and therefore R — 1/2. For a two- or three-dimensional array, it must be 
the case that at least some of the qubits are capable of interacting with more than two neighbors (or else we 
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are isomorphic to the ID case). But for each such 'additional' interaction, we must have a barrier spin to 
control it. Thus i? < 1/2 for all arrays using this architecture. If we restrict ourselves to considering simple 
regular arrays where each qubit interacts with exactly n neighboring qubits, then we readily conclude that 
the optimal value R — 2/5 corresponds to n = 3. (Notice that if we inspect Fig. 2 (b) and (c), but assign 
the opposite colour scheme to that specified in the Figure, i.e. make red— s-barrier & blue— >qubit, then these 
two geometries achieve this limit). 

It seems regrettable that moving to a higher dimensional array, with the attendant advantages in terms 
of connectivity, parallelism, etc, must apparently necessitate a lower value of R. But in fact this is only 
true if we limit ourselves to considering networks in which there is one barrier spin for each adjacent pair 
of qubits. If we can allow a barrier spin to do 'double duty', in the sense that it controls more than one 
interaction, then we can achieve higher values. For example, consider again Fig. 2(b) & (c) but now take 
the red colour as denoting a qubit, and the blue as denoting a barrier. The value of R thus increases from 
2/5 to 3/5 and we have an architecture that is more efhcient than the one-dimensional case. Similarly, the 
geometry shown in (d) would offer R = 2/3. However, the difficulty is that 'lowering' such a barrier will 
lead to a much more difficult condition in terms of the state revival: the barrier must, at some subsequent 
time, be dis-entangled from all adjacent qubits (i.e., must return to a product state) in order that it can 
be 'raised' again. We will now consider whether this condition can be met for the two simplest cases: the 
three-qubit mode, as shown in Fig. 1(b), and the four qubit node illustrated in Fig. 1(c). 

The Three Qubit Node. 

Consider a sub- network of seven spins with the topology and the initial state shown in Fig. 1 (b) , top. 
This can be seen as a 'fragment' of either the 2D pattern of Fig. 2(b) or the 3D pattern in Fig 2(c). There 
are four barrier spins in the eigenstate | t) and three qubit bearing spins, the qubits being labelled W , X 
and Y . Initially we suppose the device is in a 'passive' state. The physical distinction between the barrier 
spins and the qubit-bearing spins is simply that the two types have far off-resonance Zeeman energies, so 
that the underlying Heisenberg interaction between a barrier spin and a qubit spin is effectively of the Ising 
form. Thus the barriers will remain in their | 1) eigenstates indefinitely (or more accurately, for a finite time 
depending on the accuracy of the Ising approximation, which dependsii on the magnitude of the ratio of 
the detuning to the interaction strength \ A~ B\/J). 

Now suppose that we abruptly tune the Zeeman energy of the central barrier spin so that it is comparable 
to the Zeeman energy of the qubit-bearing spins. (The transition need not be abrupliii, but this assumption 
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FIG. 2: Illustration of the geometries that can be formed from the multi-qubit gates depicted in 
Fig. 1. Here (a) is the simple ABAB... linear array that is established in the literature. Parts (b) 
& (c) show examples of the two- and three-dimensional arrays that can be generated purely from 
a node of the kind shown in Fig. 1(b). (Fortunately, hexagonal patterns are in fact quite common 
occurrences in self-assembling nanostructures.) In contrast to the networks suggested previously, 
here the nodes of the structure are barrier spins, and therefore each such barrier separates more 
than two qubits. In (d) we see one example of the structure that could be formed from the 
four-qubit node depicted in Fig. 1(d). 

vastly simplifies the analysis). Notice that the three outer barrier spins, being still far off-resonance from 
the qubit spins, will continue to interact via an Ising form J(Jz<^z - and moreover they will remain in the 
state I X) throughout. Therefore their contribution to the Hamiltonian of the central four spins is simply a 
shift to the Zeeman energy of their neighbors. Let us label the three spins which initially bear qubits by the 
numbers 1, 2, 3, and the middle spin by the letter M . Then we can write the Hamiltonian of these spins as 

-^quadruplet — -^zeeman ~t~ -^int 
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-f^zeeman = Baf^ + {A + Jz) Cr, 



Hint = Jxvjo'fo'M + '^J'^m) + JzO'fo'M (1) 

It will prove interesting to see the consequences of planar anisotropy {Jz ^ Jxy) when we attempt to achieve 
a complete gate. A full representation of this Hamiltonian requires a matrix with 16 = 2* rows/columns. 
However it is of course block-diagonal, as illustrated in Fig. 3(a), because ifj„t commutes with Sz = 
(^M+Yli Thus we can understand the complete dynamics by considering the motion within each subspace 
separately. In the figure we introduce the notation Hx,y for the Hamiltonian within a given subspace, 
where the subscript {X,Y) indicates the space is spanned states having X 'up' spins and Y 'down' spins 
(equivalently, one could of course label by the Sz eigenvalues). We will establish the eigenstructure of these 
subspaccs, and then consider whether the 'revivals' in each space can be made to coincide. 

For the complete set of 16 spin basis states, we use a notation of the form | t) where the right 

ket contains the state of the central spin, and the left ket contains the states of the outer three spins, in 
clockwise order from the top left (although in fact the order will not matter to us). If we write the initial 
computational basis states as |000), |001)...|111), with the same ordering, then recalling that the initial 
barrier state is | t) we simply identify these initial states with the eight of the sixteen spin states thus: 

|ooo) ^ I T) 
1001) ^ I iiT)l T) 



|iii) ^ I TTT)| T) 

Firstly we note that there are two trivial subspaces associated with the states | tTT)| T) and | iii)| i) 
(though the latter is not of interest). The computational basis state |111) therefore remains an eigenstate 
with energy 3(^4 + Jz) + B + 3Jz- We introduce the symbols a = A+ Jz and S = B — a, so that this energy 
is 4a + (5 + SJz- 

Now the smallest non-trivial subspaces are those with Hamiltonians H^^i and i?i,3. As required by 
symmetry, these Hamiltonians prove to have the same structure: 

1 \ 
1 



H = K1 + 2JxY 



1 

1 1 1 X 




FIG. 3: Diagram showing the block-diagonal structure of the Hamiltonian for the three-qubit node 
(a) and the four-qubit node (b). The spin basis states are shown schematically: the notation 
I TTT)I T); for example, corresponds to the uppermost row in part (a). 

where 1 is the identity matrix. For the H^ i subspace K — K^^i = 2a + S + Jz, X = X3.1 — —{d + 2Jz)/JxY, 
and the basis runs {| jtT)l T)i I TiT)l T)7 I TTi)l T)i I TTT)I i)}- Notice that the fourth basis state is not a 
possible initial computational basis state. For the complimentary i?i.3 subspace we have K1.3 = —2a — S + Jz 
and Xi,3 - (<5 - 2Jz)/JxY, the basis running {| tii)| i), | iTi)l I), I iiT)l i), I ?)}• Notice that this 

time only the fourth basis state is a possible initial computational basis state. 

The eigenstates/ values in the two spaces of course have the same form. We may write the (unnormalised) 
eigenvectors compactly as: 
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Vectors \V±) have eigenvalues E = V± — ^ {-^ ^ (12 + X'^)^/'^) , while the paired vectors are degenerate 
with value E = 0. Here X is of course either ^3^1 or X^^^ for the two respective subspaces. From these 
eigenvalues, the full energies are obtained as K + 2JxyE. 

We are now in a position to examine how certain qubit basis states evolve. Consider the basis state 
|000). This corresponds to | t) and is therefore in the subspace of -ffi,3. 



1 000) 



= N{V+\V.}-V.\V+)) 






At time t after the Zeeman shift event, we will have: 

|000) ^ N' exp{-iT{-2a - S + Jz + 2JxyV-))(v+\V-) - V- exp{-iT2JxY{V+ - V-))\V+)^ 

where we have used r = t/U. Consider the phase difference between \V-) and |V+). If this is anything other 
than the original —1, then the state will some finite projection onto {| tii)| Di I iTi)| i)) I iiT)| i)}- Since 
these states are not computational basis states, we immediately have our revival condition: tr must be such 
that exp{-iT2JxY{V+ - V-)) = +1. This yields tr = rmr{12JxY + {S - 2Jz)'^)~^^'^ for an integer m, as 
our first revival time criterion. At such times we can abruptly shift the Zeeman energy back to its large 
detuning, and thus decouple the barrier spin. If we do so, the net effect will be simply to introduce a phase: 
|000) exp(i(/))|000) where we find ^ = TR{2a + 5 — Jz — 2JxyV-). Note that this phase in fact differs 
from the passive phase <?f>passive = t(2o — 5 + 3Jz) that would be accumulated over the same period by a 
state 1 000) if there had been no Jxy interaction. This phase difference must be remembered when we come 
to describe the process as an effective operation within the space of qubit basis states. 

Following analogous argmnents we can examine the evolution for the basis states {|011), |101), |110)} 
which are governed by H^^i. Now the revival condition is that any state within the space spanned by this 
triplet must return to it; or equivalently, there must be no projection onto the fourth, non-computational 
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basis state | tTT)| i)- In general a state \s), which at t = is orthogonal to | TTT)| i)) will develop as: 

|s) ^ A/-exp(i^i)((a|asymi) + /3|asym2)) +exp(i</)2)(y+|y_) - VI exp(i</)3)|l/+))) 

We find that our condition of zero projection onto | tTT)l i) corresponds to the constraint exp{i(j)^) = +1, 
which implies t = n7r(12j|^y + (6 + 2Jz)'^)^^^^ ■ This revival condition must be met simultaneously with 
the former condition, and indeed with the further conditions we are about to derive. 
The ^2,2 subspace remains to be characterized. We find: 



H2,2 = {6-Jz)l + 2JxY 
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Here x = -6/Jxy and our basis is {| Tii)| T), | iTi)| T), I iiT)| T), I TTi)| i), I TiT)| i), I iTT)| i)}- The 
eigenstates include two degenerate pairs: 
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with eigenvalues 14± = | (a; ± (4 + a;^)^/^). The remaining, non-degenerate eigenvectors are 
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V ^16- j 



with eigenvalues V\%± = 5 (a; ± (16 + a;^)^/^). These eigenvalues E = {14±, 146±} of course correspond to 
total energies 6 — Jz + 2JxyE. 
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Having thus determined the eigenstructure of this subspace, we can again consider the evolution of the 
corresponding computational basis states: {|100),|010),|001)}. We find that any state initially within the 
space spanned by this triplet will develop as: 

exp(i</)i) |^a(f4+ HF4-a) + /3| ^4-6)) - exp(i(/)2)y4- (a| Vi+a) + /3|F4+6))) 

+Bexp(i03)(vi6+|Vi6-) + exp(i04)Vi6-|Vi6+6)) j 

(where A, B, a, (3 arc constants depending on the initial state). 

In order for this state to be within the original space we must meet two phase conditions: exp(i02) = +1 
and exp(i04) = +1. These correspond to two constraints on the revival time: tr = P'k{AJ\y + and 
tr = g7r(16J|-y + Now we can collect the four conditions on tr and look for a common solution: 

' m(12j2y + (5-2Jz)2)-i/2 
n(12j2^ + (5 + 2Jz)2)-V2 

g(16Ji^ + 52)-i/2 

The Limit of a Pure XY Interaction (Forster Type Interaction). 

Consider the limit where Jz = 0. Then the first two constraints are equivalent, and we can rewrite the 
conditions as 

n(12 + rf2)-i/2 

-TR= \ p(4 + d2)-l/2 

g(16 + rf2)-i/2 

where we have taken Jxy to the right and introduced d = 5/Jxy- In general, there is no value of d for 
which one can find integers n, m, q to satisfy the conditions. We can find values of d such that two of the 
three roots ()~^/^ are compatible (i.e. in a rational ratio) but the third will not be. For example, = 
makes the second two compatible via p — 2q, but the first then differs by an irrational factor \/3. Thus we 
conclude that our gate cannot operate in the strict limit Jz = 0. 

The Limit of a Pure Isotropic Interaction. 

Consider instead the case where Jz = Jxy = J- In this case there is a value of 5 for which all four 
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conditions can be met: namely the trivial value (5 = 0. Then: 

m/(4J) 
n/{AJ) 
p/{2J) 

and the first revival will be — 7r/(2 J) implying m — n — q — 2 and p — 1. Given this value of r we can 
explicitly write out the transformation in the computational basis. We find: 



U = 



-10 

-1/3 2/3 2/3 

2/3 -1/3 2/3 

2/3 2/3 -1/3 





















(3) 



1/3 -2/3 -2/3 

-2/3 1/3 -2/3 

-2/3 -2/3 1/3 

^ 1 y 

The basis order here is {|000), |001), |010), |100), |110), |101), |011>, |111>}. It is interesting to note that 
= 1, thus the second revival at 2rfl simply yields the indentity: a complete return to the initial state. 
However U itself is an entangling operator and one can generate more familiar operations by combining it 
with single qubit operations. For example, one caniS obtain: 



S = 



'^lOOOOOOO^ 
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p = 



-1 000000 0^ 
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000000 -1 



where is basis order here is the same as for U. The operator S causes |010) to be exchanged with |100) 
and 1 101) with |011), thus the operation is simply a SWAP between two of the qubits (the first two in our 
notation). In terms of the arrangement in Fig. 1(b), this would correspond to swapping W and X. Of course, 
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by suitably varying the composition of this operation, we can exchange any two selected qubits among the 
three adjacent to the barrier node. Thus these operations allows us to freely permute the arrangement of 
qubits within a network such as Fig. 2(b) or 2(c). 

The operator P simply associates a phase of —1 with the states |000) and |111). This is a kind of 
extension of the familiar two-qubit gate "control-cr^" . It is interesting to note that while the two-qubit 
operation of associating —1 with |00) and |11) can be composed from single qubit gates, this three-qubit 
analogue cannot. If, for example, we were to ensure that one of the qubits in the triplet was in state |0) 
(as we could simply by swapping a known into the triplet, using S) then we can create a control-NOT 
between the other two qubits simply via XPX where X = -^{ax +<^z) is the usual single-qubit Hadamard 
operation on the target qubit. 

Thus S and P, together with single-qubit gates, constitute a universal set of gates for QC. We conclude 
that the form of barrier gate illustrated in Fig. 1(b) does indeed suffice for gating interactions. 

General Jxy ^ Jz Models. 

We have seen that the limit of the pure XY interaction {Jz = 0) is not suitable for this type of gate, 
but the limit of an isotropic Heisenberg interaction [Jz = Jxy) is suitable. Away from these limits we have 
(with d = 6/ Jxy and f2 = Jz/ Jxy) 

' m(12+(2f7-rf)2)-i/2 
n(12 + (217 + d)2)-i/2 
p(4 + d2)-i/2 
g(16 + d2)-i/2 

In general, for an arbitrarily chosen f2 there is no solution. Are there any f2 values, other than the previously 
established = 1, that will suffice? The answer is yes. For example, taking d = we find 

n(12 + 4r!2)-i/2 
p/2 
g/4 

and we see that the condition that there be a solution is simply that (12 + 4r2^)^^/^ be a rational number. 
Thus there are an infinity of suitable f2 values, and indeed for any chosen there will be an acceptable 
arbitrarily close to it. Of course, to fully assess the practicality of using this scheme for such a value, 
one would need to generate the matrix corresponding to Eqn.[3] and establish that it is entangling, find 
constructions for useful operators, etc. 
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We have seen that the three-qubit node can be used for QC, at least for a range of interactions that 
includes the important isotropic form. The natural question then becomes, can we extend the idea further 
to a four-qubit node? 



The Four-Qubit Node. 

The analysis follows exactly as the threc-qubit node, however with different conclusions. Wc consider a 
sub-network of nine spins with the topology and the initial state shown in Fig. 1(c), top. The full analysis is 
given in the Appendix. There we argue that, as with the three-qubit node, there is generally no advantage 
in considering cases where 6^0. With this simplification, the constraints for a revival are, 



JxY 
rR 



m(8 + 02)-i 

p(16 + 9r22)-i (4) 

where again = Jz/Jxy- 

In general, one can choose an to bring two of the three roots into a rational ratio - but the third will 
not be. There is only one apparent special case where all three can be satisfied, namely the = 1 limit 
which corresponds to the important isotropic Heisenberg interaction. One might think, therefore, that at 
least in this one case the four-qubit node may function. However, inspecting the resulting conditions: 

JxY 



IT 

Jxy 



-tr = m/3 = p/5 = s/5 



we find that the only solution is to set ^^^tr equal to an integer, so the first revival is with m = 3, p = s = 5. 
When we enter this value of tr into the explicit expressions for state evolution, we find that the net effect 
of the entire process is simply the identity matrix. Thus, whereas for the three-qubit node we found that 
there was a useful revival at r = tt/ (2 J) prior to the second, identity revival at r = n/J, in the present case 
of the four-qubit node we find that there exists only the identity revival at r = tt/J. Therefore we must 
conclude that it is not possible to implement a useful four-qubit gate in this way. 

To summarize, we have examined two forms of multi-qubit gate based on a barrier scheme for QC which 
permits interactions to remain 'always-on'. For the three-qubit node, we found that the isotropic Heisenberg 
interaction indeed yields an entangling gate operation on a time scale of n/{2J), and confirmed that this 
operation can be used to construct more familiar SWAP and Control-NOT operations. Moreover, it seems 
possible that most, if not all, anisotropic Heisenberg interactions can be supported to a sufficient level of 
accuracy for QC. However, for the four-qubit node, we found that no useful operation can be performed. 
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We have not ruled out the possibihty that manipulating the qubit-bearing spins during the gate operation 
(by shifting Zeeman energies or by rapid flipping, for example) might generate a useful net operation, but 
this is beyond the scope of the basic barrier scheme. Possibly this is a direction for future work, as is the 
problem of finding optimal constructions for canonical gates via the thrcc-qubit node. 

We conclude that the network topology shown in Fig. 2(d) is not feasible, but that the forms of network 
shown in Fig. 2(b) and (c) are feasible. These networks have a qubit storage efficiency that is 50% higher 
than previously proposed 2D and 3D schemes. Moreover these networks are more efficient than a one- 
dimensional array, and may therefore represent the 'best yet' embodiment for quantum computing with 
always-on interactions. 

The author wishes to acknowledge support from a Royal Society URP, and from the Foresight LINK 
project "Nanoelectronics at the Quantum Edge" . 
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Appendix: Analysis of the Four-Qubit Node. 



Consider a sub-network of nine spins with the topology and the initial state shown in Fig. 1(c), top. Now 
there are five barrier spins in the eigenstate | f) and four qubit bearing spins, the qubits being labelled W, 
X, Y and Z. As before we tune the Zeeman energy of the central barrier spin so that is comparable to the 
Zeeman energy of the qubit-bearing spins, and we assume for simplicity of analysis that the change is abrupt. 
By the same reasoning stated previously, the can neglect the outermost barriers except for the Zeeman shift 
they create. Then labelling the four spins which initially bear qubits by the numbers 1..4, and the middle 
spin by the letter M, the Hamiltonian for these five spins is exactly like Eqn.pQ except that the index of 
course runs j — 1...4. Complete representation of this Hamiltonian requires a matrix with 32 rows/columns, 
with the block-diagonal structure illustrated in Fig. 3(b). As before we will establish the eigenstructure 
of these subspaces, and then consider whether the 'revivals' in each space can be made to coincide. The 
notation is analogous to that used before: | t), etc., where the left ket contains the states of the outer 

four spins, in clockwise order from the top left. Similarly we write the initial computation basis states as 
|0000), |0001)...|1111), where the order is again clockwise from the top left, making the association exactly 
analogous to Eqn.|2]. 

The two trivial subspaces are now associated with the states | tTTT)l T) | J,iii)| J,), which therefore 
remain eigenstates with energies AJz ± (5a -|- S). 

The smallest non-trivial subspaces are now those with Hamiltonians H^ ^ and -ffi,4. Again, these Hamil- 
tonians will have the same structure due to symmetry. The form is: 

1 ^ 

1 

1 

1 

1 1 1 1 a: 



H.i ^Kl + 2JxY 



For Hij we have K ^ {3a + S + 2Jz) and X = X^s = -{S + 3Jz)/Jxy, and the basis is {| iTTT)l T), 
I TiTT>l T>, I TTiT)l T), I TTTi>l T>, I TTTT)I i)}- For i/1,4 we have K - (-3a - S + 2Jz) and X = X^s = 
{S - 3Jz)/JxY, and the basis is {| Tiii)| i), | iTii)| i), | iiTi)l i>, I iiiT>l i), I iiii>l !)}• 
We may write the (unnormalised) eigenvectors as: 
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|ai) 



V / 



|a2) 





1 

-1 

V y 



las) 



V / 



|V-> 



|V+> 



1 
1 
1 



The vectors |y±) have eigenvalues E = V± = ^ [X ± (16 + X^)^/^), while the bracketed triplet are degen- 
erate with value E = 0. Here X is of course cither X4.1 or ^1,4 for the two respective subspaces. From 
these eigenvalues, the full energies are obtained as iiT + 2JxyE. By the same reasoning employed for the 
three-qubit node, we find one revival time constraint from each of these subspaces. From H4^i (where the 
initial state will be in the space spanned by the first four spin states, and at revival must return to this 
space) : 

TR = mn{16J^Y + + ^Jzfr^ 
From -ffi,4, where the initial state is |0000) = | t) and at revival must return to this state: 

TR = n7r(16Jiy + {S- SJzfyi. 

The remaining two subspaces, if3_2 and -ff2,3i again have a common structure: 

'^OOOOOOOOll^ 
0000000110 
0000001100 
0000001001 
0000001010 
0000000101 
001110a;000 
0110010a;00 
11001000x0 
1 1 1 a; 

For H3^2 we have K = {a + 6) and x = ^3,2 = -{6 + Jz)/Jxy, and the basis is | iiTT)| T), | TiiT)| T), 
I TTii)| T)> I iTTi)| T)> I TiTi)| T)> I iTiT)| T) for the first six (the computational basis states), and 
I TTTi)| i), I TTiT)| i), I TiTT)| i), I iTTT)| i) for the remaining four. For H2,3 we have K = -{a + S) and 
X = X2,s = — Jz)/JxY, and the basis is of course the compliment of the Hs^2 basis. 



H = K1 + 2JxY 



\ 



J 
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The eigenstates include a degenerate pair {lasym^), |asym2)} with eigenvalue E = and two non- 
degenerate vectors IV24+), IV24-) with values E = V2i± = \(x± (24 + a;^)^/^): 



|asymi) 



I 1 \ 

-1 
1 

-1 






V / 



lasyms) 



' 1 ^ 



1 
1 
1 

-2 
-2 




V / 



|V24±) 



The remaining eigenstates form two degenerate triplets, 
{\Vi+a), \Vi+h), \V&+c)}, with eigenvalues E = Vs± = \{x±{S + x'^fl'^). 



|V8±a) = 



/ 1 ^ 



1 

-1 
-1 





y%± 




V 



1 

-1 
-1 

1 











\Vi±c) 



( 2 \ 
2 
2 
2 
2 
2 

V24± 
^24± 
V^24± 

V"24± j 

{|l/8-a), 11/8-6), |14-c)} and 

/ ^ 




2 

-2 



all eigenvalues corresponding of course to total energies 5 — Jz + 2JxyE. 

We can now consider the motion of states initially in the computational subspaces. First consider the i?2,3 
space: here the computational subspace is {|0001), |0010), |0100), |1000)}. By analysing the development 
of a state initially in this space, exactly as was done for the three-qubit node, we find two constraints on 
the revival time: 

mr{2AJ%Y + - Jzfy^ 
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Similarly, from the H^ 2 space containing the six computational basis states {|1100), |0110), |0011), |1001), 
|1010), |0101)}, we find that the condition for a state, initially within that space, to retmm to it is: 



TR 



nn:{SJj.y + {5 + Jzf)--^ 
mi{2AJly + {5 + Jzf)-l 



J 



XY 



-TR 



Combining these conditions with those from the smaller i/4^i & -ffi,4 spaces, we have our complete set of 
revival time conditions: 

m(8+ + 
n(8 + (r2-d)2)-^ 
p(16 + (3fi + (i)2)-i 
g(16+(3f7-rf)2)-i 
r(24+ + 
s{2A+{n-d)f)-^ 

As before, we observe that there will exist some set of integers m..s which indeed mutually satisfy these 
conditions if, and only if, all ()^^ terms are be related to one another in rational ratios. We note that there 
is no advantage in using d ^ Q since this doubles the number of conditions while only providing one variable. 
Therefore using d — Q we find 

p(16 + 9r22)-i 
s(24 + r22)-i 



Jx\ 



-TR 



as quoted in the body of the paper, Eqn.0]. 



^ NMR systems are a partial exception to this rule: there the interaction is of the diagonal Ising 
form and is weak compared to the rate at which single spins can be rotated; the combination of 
these two special factors allows the 'always on' interaction to be effectively negated by methods 
such as periodic spin flipping. 
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A one-dimensional architecture that relies on the nearest-neighbor interactions must pay a heavy 
overhead in terms of the proportion of operations that are simple SWAPs, needed merely to 
permute qubits. This is due to the average separation between a given qubit pair being oc in 
a ID array, versus A^/^ or N'^^^ in two- or three-dimensional networks. There are consequences 
for fundamental considerations such as the tolerable error rate etc. 

^3 D. Loss k D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 
R. Vrijen et al, Phys. Rev. A 62, 012306 (2000). 

15 A. Imamoglu et al, Phys. Rev. Lett. 83, 4204 (1999). 
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J. Siewert et al, J. Low. Temp. Phys. 118, 795 (2000). and references therein. 
The following argument will allow the reader to verify our claim that S and P can be 
formed. These constructions are in no sense optimal - finding the efficient compositions 
for specific gates is beyond the scope of the present paper. Let Rz{a, (3,^) be the opera- 
tor that applies a single qubit Z-rotation cos all -|- i sin aaz to the first qubit, while inde- 
pendently rotating the second and third by (3 and 7 respectively. Introduce special cases 
R± = Rz{x,-x,0) with X = ±(7r/4 arcsin(l/8)/2). Let U± = U.R±.U. Then defining 
W = R{-TT/2 + y,-TT/2-y,-TT/2).U+.R{0,0,TT/2).U- where y = arctan(-^/7/3)/2, we notice 
W = P.S i.e. it is composed of the two operators we are interested in. Further, we find that 
R{7r/4,7r/ 4, 0).W.Tyz-W.Tyz-T+ = P, where single-qubit gates Tyz = ^{<^Y + <^z) and 
T± = :^(1 i icTy) all act on the third qubit. Given P we know we must also be able to obtain 
S since S = P.W. 
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